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The new high precision data on charmonium production in proton- nucleus collisions by the E866/NuSea collab- 
oration at Fermilab allow — together with older data at lower energies — to fix a unique set of parameters for 
the standard production and absorption scenario of charmonium in a proton-nucleus reaction. In this scenario the 
cc pair is formed in an octet state, emits a gluon and continues its radial expansion in a singlet state until it has 
reached the charmonium radius. In all three phases it can interact with the nuclear environment. We find that 
the lifetime of the octet state is much shorter than acceptable on physical grounds. This challenges the physical 
reality of the first phase in the standard scenario. 
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I. INTRODUCTION 

Since the suppression of J/ip's has been proposed by 
Matsui and Satz as a signal for a quark-gluon plasma 
formation more than a decade ago Q a lot of experi- 
mental efforts have been devoted to use them for the 
quest of the formation of a quark gluon plasma created in 
ultra-relativistic heavy ion collisions. For small systems, 
where the creation of a plasma is not expected, the mea- 
sured absorption can be well explained without invoking 
a plasma. At Quark Matter 1996, the first results on the 
heavy system Pb-Pb have been reported [|j by the NA50 
collaboration and it has been argued that they show an 
abnormal suppression (as compared to the extrapolation 
from small systems) which may be interpreted as a hint 
for the formation of a quark gluon plasma. This triggered 
a lively and still ongoing debate about the conclusions of 
this observation. 

Before one tries to understand the J/ip suppression 
seen in heavy ion collisions it is necessary to understand 
that obtained inpA collisions where the reaction is much 
more under control. There one expects that the shape 
and the density of the nucleus does not change during the 
reaction. Recently new data have been advanced by the 
Fermilab group E866/NuSea j^, where for the first time 
differences between the absorption of J/tp^s and tp'''s have 
been observed. This fact as well as the high precision of 
these data allow to fix the parameters of the standard 
scenario 

Here, we present a model for the charmonium ab- 
sorption in proton-nucleus collisions which is based on 
the standard scenario described in Section II. We have 
limited the ingredients to the minimum that allows the 



description of the Fermilab data in the whole x range. 
Several initial state effects such as gluon shadowing [p|, 
energy loss 0,1), or quantum coherence [D could af- 
fect somehow charmonia production at Fermilab energies. 
However, since Drell-Yan production does not show a sig- 
nificant nuclear dependence, we have assumed that those 
effects are relatively less important and their considera- 
tion would imply an unnecessary complication. We have 
checked that inclusion of gluon shadowing does not alter 
the conclusions one can draw from our study, but changes 
somewhat the present set of effective parameters. Tak- 
ing the recent results of the E866/NuSea collaboration 
at Fermilab ^ , we determine a set of parameters for the 
standard scenario. The same set of parameters allows the 
description of almost all other available data on charmo- 
nium suppression inpA collisions, including the majority 
of NA3 data points (lO|, those of NA38 ||ll| and those of 
E772 ||l^ . We show that the functional form of the model 
allows to describe the whole body of data and that its pa- 
rameters are rather precisely fixed by them. We show as 
well that the X2 scaling predicted by the model is com- 
patible with high-energy data. We do not claim that our 
result excludes other reaction mechanisms (as we will see 
there is good evidence for that on physical grounds) but 
stress that on the basis of the presently available data 
one cannot distinguish experimentally between them. 



II. MODEL FOR CHARMONIUM ABSORPTION 
IN PROTON-NUCLEUS REACTIONS 

InpA collisions the production of cc states is supposed 
to be proportional to the atomic mass A of the target, as 
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it is the case for the Drell-Yan process. Thus, the charmo- 
nium states are produced anywhere in the nucleus with 
a probabihty directly proportional to the local nucleon 
density. We consider here J/tp^s, ■0"s and x's. Feeding is 
important: only 50% of the observed J/i/j's are primarily 
J/'^'s, the other 50% come from radiative decays of x's 
(43%) and V"s (7%) §. The spin of the x state may 
play a role in the x ^ cross section , but we neglect 
this effect here. 

The cc pair is not created in a fully formed state, with 
its final radius. We assume that the transverse distance 
between the quarks increases linearly with time. In the 
cc rest frame we write 



ro + Vc 
n 



if rce(r) < ri 
otherwise. 



(1) 



where i stands for J/ip, ip' and x- The final radii of 
the different states are determined by microscopic calcu- 
lations. The limits chosen are 0.43 fm for J/V', 0.87 fm 
for ip' and 0.67 fm for x il4| . Wce is the sum of the c and c 
velocities i.e. Vcs = 2uc. fee and the initial radius tq are 
taken as free parameters. An alternative possibility for 
the cc expansion, r (x \/t, has been proposed in |15|. We 
have checked that no definite conclusion may be drawn 
about which one is actually at work. The main point is 
that the expansion is important in the kinematical region 
where one can distinguish J/i/; from V''. 

The cc pair is formed by two-gluon fusion in an octet 
state (cc)8 which later on emits a gluon and neutralizes its 
color. We assume that the emission occurs after a color 
neutralization time rg^i in the rest frame of the cc pair, 
independently of its velocity, rg^i is the third parameter 
of the model. In the nucleus rest frame, the lifetime of the 
octet state increases with xp, due to Lorentz dilatation. 
We will see that this fact is responsible for the increase 
of absorption with Xi?, as seen in data. This would not 
be the case if the octet-singlet transition occurred after 
a constant time in the nucleus rest frame (a natural hy- 
pothesis if neutralization was caused by interactions with 
the surrounding nucleons). According to the octet 

lifetime should be of the order of 0.3 fm. We will see 
that the new data allow for a precise determination of 
the color neutralization time rg^i. 

Some authors have argued that the production can 
partly proceed via a direct color sing let state Al- 
though a priori distinct from color evaporation they have 
shown that it is not presently possible to clearly disen- 
tangle those two possibilities. Since on the one hand the 
xp dependence of the color singlet production fraction is 
unclear and on the other hand the fit does not improve if 
one incorporates this effect we assume that cc pairs are 
always created in a color octet state. 

The singlet-nucleon absorption cross section cr^^jgj^Tv 
depends on both the transverse radius r^c of the cc pair 
and its energy. In QCD the total cross section for a com- 
pact singlet state should be proportional to the square of 
its radius and we write oi^^^-^^j^ — ct^n{s) ■ [tcc/t^Y. 



Next, we need the energy dependence of the tpN cross sec- 
tion in the kinematical region corresponding to available 
data, that is in the region around 10 GeV and above. 
In this range it may be related to that for J/?/' photo- 
production and to the small x gluon distribution M] 
leading to cr^Ar(s) = (Ti • (y^/10 GeV)°-^ where ^/s\s 
the center of mass energy of the ccA^ system. Thus, the 
time dependent cross sections for singlet states is given 
by 



(^(cc)in{t) = CTl 



10 GeV 



(2) 



where the radius r^g is given by and cti is the cross 
section for a fully formed J/ip with incident energy ^/s — 
10 GeV. CTi is the fourth parameter of our model and 
should be of the order of a few mb 1 1^ . 

The present model gives a unified treatment for J/ip, 
tp' and X since the only difference in the cross sections 
comes from the final radii r^. Gonsequently, one expects 
for small xp where the J/V' is formed before leaving the 
nucleus a difference in the absorption between J/tp, ip' 
and Xj as the V'' and the x radii are bigger than that of 
the J/tp- Indeed the E866/NuSea group found a different 
absorption tor xp < 0.2 (cf. Figs.]^ and||). At high xp, 
the different states will behave similarly because the cc 
leaves the nucleus before it has reached the J/tp final size, 
due to Lorentz dilatation. Strictly speaking this is true 
if one neglects the mass difference between charmonia. 
For given beam energy and charmonium velocity xp is 
proportional to the charmonium mass implying, e.g., a 
systematic 20% shift of tp' with respect to J/tp. This 
is somewhat reduced by the fact that J/tp comes partly 
from radiative decay of x and tp' and further by the fact 
that what one needs in order to compute the cc velocity at 
a given xp is the intermediate cc mass. In lack of a good 
description of this mass effect and noticing that present 
tp' data are not precise enough to be analyzed at this level 
of accuracy, we assumed that the intermediate cc and all 
charmonia have the same mass rricc = (m^, + m^i )/2 = 
3.4 GeV. 

The interaction between the octet state and nucleons 
is mainly responsible for the increase of absorption with 
increasing x p observed at Fermilab. The large x p region 
is dominated by cc pairs which have left the nucleus in the 
octet state. The absorption cross section between octet 
states and nucleons (J[cc)aN rnay be of the order of 20 
mb 1^, that is much larger than the singlet cross section. 
The usually quoted value of 6 mb in hadroproduction 
would be interpreted in the present scenario as an average 
between singlet and octet cross sections. Since the energy 
dependence is related to gluon dynamics, we assume that 
it is the same for CT(cc)sAf and <J(^cc)iN- The octet cross 
section is taken to be independent of the radius of the 
pair (cc)8: 



10 GeV 



0.4 



(3) 
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CTg is the fifth and last parameter of the model. 

If one assumes singlet states only, the A to p produc- 
tion ratio approaches 1 at xp = 1 because the large 
gamma factor between the cc pair and the laboratory 
frame does not allow an expansion of the cc system while 
traveling through the nucleus. Thus the cross section (H) 
is small and so is the absorption. 



III. RESULTS 

Quantitative results were obtained with a Monte-Carlo 
simulation. The computation is done by first generating 
3 A random numbers, the first 3{A — 1) describing the 
positions ot A — 1 static nuclcons inside the nucleus and 
the 3 last localizing the creation point of the charmonium 
inside the nucleus. From this creation at the proper time 
T = 0, we follow the charmonium on its way through the 
nucleus in the z direction, z(t) = zq + (3^t. The charmo- 
nium velocity (3 in the nucleus rest frame is determined 
hy xp (see Section IV). The charmonium neutralizes its 
color at time rg^i and expands as expressed by (^. If 
the transverse distance between the charmonium and a 
nucleon becomes smaller than ■y/(T(r)/7r we assume that 
the charmonium is absorbed. 

We made a fit to the J suppression in tungsten nor- 
malized to beryllium, measured by the E866/NuSea col- 
laboration 1^, with the five parameters described in the 
previous section. We obtain the following parameters: 
ro = 0.15 fm, Vcs = 1-85, cti = 2.1 mb, rg^i = 0.02 fm 
and CTs — 22.3 mb. This set of parameters is used to cal- 
culate the results for all other energies and targets. The 
influence of a change of these parameters on the results 
is discussed in Section VI. 

In Fig. |l| we compare our calculation with the re- 
cent E866/NuSea results of proton- nucleus reactions at 
800 GeV. In this experiment the absorption ratio 



R{A/Bc) = 



9g(A) 
ylcr(Be) 



(4) 



for W and Fe has been measured. First of all we see 
that below xp = 0.2 the absorption of J/ijj and tp' is dif- 
ferent. It is this fact which allows to determine the pa- 
rameters of our model precisely. We see as well that the 
absorption has a minimum around xp — QA. As we will 
discuss later in detail the increase at large a;^ is caused 
by octet state-nucleon collisions whereas that at negative 
X F comes from collisions of the fully formed charmonium 
states with nucleons, thus explaining the difference be- 
tween J and absorption. The above parameter set 
describes data quite well in the whole measured x p range 
and for the two different ratios (notice that the fit was 
performed with the W to Be production ratio only). 

In Fig. H we show the same result in the alternative 
representation 



a = 1 



lnfl(W/Be) 
In 184/9 ■ 
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FIG. 1. Ratio _R(yl/Be) = {9 /A) ■ cr(A)/cr(Be) versus xp of 
produced J/ip's (black circles) and tp''s (white circles) for W 
(top) and Fe (bottom). The lines are the results of our model. 
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A thorough study of this quantity |Q has shown that 
E866 /NuSea results are compatible with earlier less pre- 
cise results at 800 GeV and for various nuclei We 
present them here for completeness though in contradis- 
tinction to the E866/NuSea data they are not corrected 
for pt acceptance [g^l ■ 
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FIG. 2. a = 1 + ln(J?(W/Be))/ln(184/9) versus xf from 
E866/NuSea for 800 GeV/c protons as compared to the re- 
sults of the calculation. E772 data are also shown for com- 
parison. 

The NA38 collaboration measured the J/ ip production 
for several targets (C, Al, Cu and W) at a beam energy 
Ep ~ 450 GeV The center of mass rapidity range 

is [—0.4,0.6] corresponding to —0.1 < xp < 0.15. The 
comparison with our model is displayed in Fig. ^ with 
xp taken as 0.05. 

Finally, we compare our model with NA3 data [ p^ . 
Measurements were done with proton as well as pion 
beams on a proton and on a platinum target, at 200 GeV. 
Here, one can compare data with the model for < xp < 
0.6. One observes in Fig. ^ a fair agreement except for 
the two largest xp values. 

For the largest xp value the deviation is large and the 
trend shown by large x p NA3 data points cannot be de- 
scribed in our model. This may be the signal that some 
energy loss of the incoming gluon manifests itself since a 
loss as given in Rcf. |Q would indeed have more influence 
at smaller energies and larger x p . We found that the sup- 
pression seen for the large xp NAB data points may be 
explained with an energy- independent 1 GeV/fm energy 
loss of the incoming gluons in addition to the nuclear sup- 
pression considered so far. We have checked that such a 
rate gives a negligible effect at 800 GeV and that the cor- 
responding rate for quarks is compatible with Drell-Yan 
measurements at 280 GeV Ell and 800 GeV 
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FIG. 3. J/tp absorption for various nuclei in p(450 GeV)+A 
reactions from NA38 as compared to the calculation. 



There are nevertheless some weak points in this expla- 
nation: 



the energy dependence 
clear 



of energy loss is not 



• the energy loss also affects the ratio at low xp for 
200 GeV which becomes more suppressed than seen 
in the data; 

• the pion beam NA3 data do not show a strong sup- 
pression of J/tp^s at large xp; 

which all together prevent us from drawing a firm con- 
clusion on the relevance of energy loss. 

Is there Scaling? Once a good reproduction of the raw 
data has been achieved, it is natural to ask whether a 
more unified picture is possible, exploiting scaling effects. 
In the present scenario, the absorption of charmonium 
states is entirely determined by the 7 factor of the cc 
pair in the nucleus frame. Therefore, the model implicitly 
contains a 7 scaling. 

The value of 7 is directly related to the momentum 
fraction X2 in the target nucleon of the gluon which pro- 
duces the cc pair. From 7 = Ecs/nicc, Ecc = xiEp and 



2xiX2'mpEp we find 



7 = 



2x2mp 



(6) 



With rricc fixed as explained in Section II, the scaling 
with 7 thus yields a scaling with X2- 

A comparison of the E772 (which are in agreement 
with the more precise E866/NuSea data) with the NA3 
data led to various conclusions in the literature. In [E3| a 
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FIG. 4. Ratio of J/?/)'s produced in p(200 GeV)+Pt and 
p(200 GeV)+p versus xp from NA3 as compared to the cal- 
culation. 

good evidence was found for X2 scaling using pion beam 
data, whereas in scaling was found in xp and not in 
X2. Since NA3 and E866/NuSea suppression ratios are 
consistent for X2 > 0.06, all conclusions depend on the 
small X2, i.e. large xf, NA3 pA data. 

IV. THE SIGNIFICANCE OF THE DIFFERENT 
KINEMATICAL REGIONS 

We now discuss in more detail which of the above- 
mentioned processes affect the suppression at a given xp 
value. 

Since rg^i is smaller than the formation time, one can 
distinguish three cc states: (a) a cc still in a color octet 
state, (b) a cc already in the singlet state but still ex- 
panding and finally (c) a fully formed charmonium state. 
Fig. ^ shows in which state the surviving cc pairs leave 
the nucleus as a function of xp in the reaction p(800 
GeV)-l-W. We see that at large xp the octet fraction is 
large whereas at negative xp values almost all cc pairs 
have lost their color at this stage. 

A complementary information is provided in Fig. |^. 
Here we display in which state the non surviving cc pairs 
are absorbed, again as a function ofxp. This figure shows 
which kinematical region is sensitive to the different cross 
sections and formation times. At large xp values the ab- 
sorption rate is sensitive to erg only. At intermediate xp 
{xp G [0,0.2]), one encounters not fully formed singlet cc 
states. The corresponding singlet cross section is weak, so 
is the associated suppression. Therefore the absorption 
has a minimum in this region. For xp < 0, the absorp- 
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FIG. 5. Relative fraction of the different states, (a), 
(b) or (c), of the cc pairs when leaving the nucleus in 
p(800 GeVHW. 

tion is mainly governed by states in expansion and fully 
formed states. Almost 50% of J ftp's are absorbed in a 
fully formed state at xp w —0.1. Therefore, differences 
between the charmonium states can be observed in this 
kinematical region. 

V. SENSITIVITY TO THE FIT PARAMETERS 

Our fit yielded a minimum at Tg^i = 0.02 fm. In or- 
der to quantify this observation we performed several fits 
in which rg^i was fixed and the other parameters were 
allowed to vary freely. In Fig. |^ the x^/ndf for differ- 
ent values of rg^i are presented. We see a rather sharp 
minimum around Tg^i = 0.02 fm. Thus we can conclude 
that the data determine the octet lifetime precisely. 

We show in Fig. ^ and || how precisely data determine 
the other parameters for rg^i = 0.02 fm. As discussed 
in the former section we can almost separate two regions: 
Xp > 0.4 where the physics is determined by the proper- 
ties of the octet state and xp < where the singlet state 
dominates. Fig. ^ shows the dependence of our results on 
(Tg. If (Tg decreases the absorption becomes weaker. This 
is a dramatic effect for large xp, which in consequence 
determines this cross section quite precisely. 

The variation of our results for J/ tjj and ip' production 
as a function of the three parameters which describe the 
interaction of the singlet state is shown in Fig. ||. Here xp 
is fixed to —0.1. We see that all these parameters are not 
very precisely determined. A change of the parameters by 
50% changes the results by about 5%. Data at lower xp - 
for both ratios R{W/Be) and R{Fe/ Be)- would therefore 
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FIG. 6. Relative fraction of the different states of the 
cc pairs at the time of absorption by a nucleon in 
p(800 GeV)+W. 
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FIG. 7. xVndf for different rs^i. The number of degrees 
of freedom is 35. 
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FIG. 8. Influence of ag on charmonium absorption in W 

and Be. The ratio of the calculation and the E866/NuSca 
data is shown for J/ip at xf = 0.12 (left) and xf = 0.82 
(right). 



be very helpful to fix them more precisely. 
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FIG. 9. Influence of ro and Vcc on charmonium ab- 
sorption in W and Bo. The ratio of the calculation and the 
E866/NuSea data is shown for J /ip (left) and tp' (right) at 
xf = -0.065. 

The initial radius ro and the quark-antiquark relative 
velocity Vcc were given by the fit as ro = 0.15 fm and 
Vcc = 1-85. Accordingly, the time needed to form & J/ili 
is 



Tf = 



r y, - rp 

Vcc 



0.15 fm. 



(7) 



For a longer formation time the ip'^s would behave as the 
J/'0's even at negative xp values, in contradiction to the 
experimental results. A shorter formation time fails to 
explain the rise of a between —0.1 < xp < 0.1. Calcu- 
lations based on realistic potentials for quark-antiquark 



6 



bound states give a somewhat higher value (0.3 fm) for 
the time of formation fE4i. 



VI. COLOR OCTET LIFETIME 

The fitted value Tg^i is very small. In order to un- 
derstand the origin of such a short time in the present 
approach we concentrate on the large x f Fermilab data. 
Here the suppression is important and absorption in 
the singlet channel with a cross section of at most a 
few mb plays a marginal role. Then the essential as- 
pects are the octet cross section and lifetime, and more 
precisely the comparison between the mean free path 
^8 = (po o'(cc)aN)~^ and the length traversed in nuclear 
matter while in octet state. The latter is bounded, on 
the one hand, by 7 Tg^i (the velocity is almost 1 in the 
region we consider) and, on the other hand, by twice 
the nuclear radius, 2Ra- There are then two extreme 
regimes: 

• 7x8^1 > 2Ra- In this regime the suppression is 
dictated by Ra/Is and the xp dependence comes only 
from the energy dependence of a(^cc)aN (see Section II). 
Since the ccN center of mass energy is 2xiEpmp and 
xi w xp at large xp^ we expect a slow dependence: 
Ra/Is (xi?)"'^. This behavior is too weak to describe 
the decrease seen at large xp. 

• 2Ra 3> jTs~,i. Here the suppression is set by 
7x8^1/^8- With 7 = xiEp/nicc, the xp dependence is 
now jTs^i/ls oc {xp)^'"^. This comes close to the ob- 
served value. The intermediate regime is an interplay of 
both scales. 

In turn the relation between the two scales may help 
to fix a maximum octet lifetime by a qualitative look at 
data. We see in Fig. 1 that the xp dependence seen in 
data does not show any clear leveling off except maybe 
in the region xp > 0.6 in both W to Be and Fe to Be 
ratios. Taking this 0;^ as a conservative lower limit for 
the transition we deduce an upper bound for t^^i: 

Tg^i < 2Rfc/i{xp = 0.6) = 0.06 fm. 

Notice that using the same reasoning with the fitted 
value of T8^i one may estimate the value of 7 or X2 at 
which one sees saturation for a given A. For W we find 



^2 I transition 



rriccTs- 



2rapR\N 



We may also go one step further and estimate the 
regime at work from the actual xp dependence seen in 
data. The variation of the number of octet states in a 
slice of matter may be written as dN = —Ndz/ls, i.e. 
N{z) — N{zo) exp[— (z — zo)/ls]- Combining this indica- 
tive exponential behavior with the xp dependence ex- 
plained above for the relevant length scale and ^8, the 
octet state suppression may be parameterized as 

5(xj.) = 5^""/""°^", 



where n depends on the regime at work. The value of n 
may be determined from data by performing the follow- 
ing evaluation 

S'jxp^) 

xPq , 

bo In bo 

where S' is the derivative of S with respect to a;i?. At 
xpg — 0.6 one finds n ^ 1 which ensures that the xp 
dependence cannot be due to the slowly varying cross 
section. 

The value of Tg^i necessary to reproduce data is very 
short compared to several theoretical estimations [ p5[ . 
Such a small value has been advocated by Wong |17| ] 
referring to an hybrid scenario proposed by Kharzeev and 
Satz in order to cure these conceptual problems. 

The short lifetime may question, however, whether the 
proposed scenario is realistic at all. The corresponding 
minimal energy of the gluon emitted for color neutral- 
ization is of the order of tens of GeV, which makes no 
sense. Therefore it may very well be that the octet state 
is just a parameterization of a much more complicated 
process. The suppression at large xp may be attributed 
to a dispersion of energy in transverse direction. Such an 
effect is expected if one assumes that the J ftp's are not 
produced directly by gluon fusion but by fragmentation 
of a color string formed between the nucleons, taking into 
account the string interactions with surrounding nucle- 
ons (the probability to have these final state interactions 
depending on the length of the path the string travels 
inside the nucleus). This alternative is currently under a 
more detailed study. 



VII. CONCLUSION 

All presently available data including recent results 
from the E866/NuSea collaboration on J ftp production 
in p A collisions can be well described in the standard 
color neutralization and state expansion scenario with a 
common set of parameters. This is the result of a Monte- 
Carlo based model in which this scenario is confronted in 
detail with data. 

It is questionable, however, whether this model is, as 
far as the color octet state is concerned, more than a well 
chosen set of fit parameters. The octet lifetime which 
can be precisely determined by the new data is too short 
for being understandable in terms of a physical process. 
Therefore it is probable that the physics observed at large 
Xp is more complex than assumed up to now. 

This parameterization of the production and absorp- 
tion of charmonium in a hadronic environment could be 
extended to nucleus-nucleus collisions. 
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